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RELATIONS BETWEEN INVASION PERCOLATION AND 
CRITICAL PERCOLATION IN TWO DIMENSIONS 

By Michael Damron 1 , Artem Sapozhnikov 2 and 
Balint Vagvolgyi 3 

Courant Institute, EURANDOM and Vrije Universiteit Amsterdam 

We study invasion percolation in two dimensions. We compare 
connectivity properties of the origin's invaded region to those of (a) 
the critical percolation cluster of the origin and (b) the incipient 
infinite cluster. To exhibit similarities, we show that for any k > 1, the 
fc-point function of the first so-called pond has the same asymptotic 
behavior as the probability that k points are in the critical cluster 
of the origin. More prominent, though, are the differences. We show 
that there are infinitely many ponds that contain many large disjoint 
p c -open clusters. Further, for k > 1, we compute the exact decay rate 
of the distribution of the radius of the fcth pond and see that it differs 
from that of the radius of the critical cluster of the origin. We finish 
by showing that the invasion percolation measure and the incipient 
infinite cluster measure are mutually singular. 



1. Introduction. Self-organized criticality has become a subject of great 
interest in recent years. Although there is no general definition for it, we 
can say that a system or model has this property if the definition of the 
model requires no parameter, yet some characteristics of the model resem- 
ble those at criticality of a parametric model with a phase transition. One 
such model is invasion percolation, a stochastic growth model that mirrors 
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aspects of the critical Bernoulli percolation picture without tuning any pa- 
rameter. The invasion model was introduced independently by two groups 
([2] and [11]), who studied it numerically. The first mathematically rigor- 
ous study of invasion percolation appeared in [4]. Connections between the 
invasion cluster and critical Bernoulli percolation have been established in, 
for instance, [4, 6, 18, 21] and [22], using both heuristics and rigorous argu- 
ments. These results indicated so many parallels between the invaded region 
and the incipient infinite cluster that a question naturally arose: to what 
extent are these objects similar? This question was studied on the regular 
tree in [1]. It was shown that, although the invaded region and the incipient 
infinite cluster are locally similar, globally, they differ significantly In this 
paper, we prove local similarities between critical Bernoulli clusters and cer- 
tain invaded clusters (the ponds) in the plane. We also show that, globally, 
the invaded region and the incipient infinite cluster are essentially different. 

In the remainder of this section, we define the invasion percolation model 
and, using results of [4], we introduce the ponds of the invasion. We then 
review results concerning relations between invasion percolation and critical 
Bernoulli percolation. Finally, we state the main results of the paper. 

1.1. The model. For simplicity, we restrict ourselves here to the square 
lattice. Invasion percolation can be similarly defined on other two-dimensional 
lattices and the results of this paper still hold for lattices which are invariant 
under reflection in one of the coordinate axes and under rotation about the 
origin by some angle in (0, tt). In particular, this includes the triangular and 
honeycomb lattices. 

Although our results concern invasion in the plane, we give the definition 
of invasion percolation for Consider the hypercubic lattice Z rf with its set 
of nearest neighbor bonds E rf . We denote edges by their endpoints, that is, 
we write e = (x, y) if the two endpoints of e are x and y. Letting G = (V, E) 
be an arbitrary subgraph of (Z d ,E d ), we define the outer edge boundary AG 
of G as follows: 

AG={e= (x,y) € E d : e ^ £(G), but x€G or j/GG}. 

The first step is to assign independent random variables, uniformly dis- 
tributed in [0,1], to each bond e € E d . We denote these variables by r e . 
Using them, we recursively define an increasing sequence Go, Gi, G2, ■ ■ ■ of 
connected subgraphs of the lattice. Go only contains the origin, with no 
edges. Once Gj = (Vi,Ei) is defined, we select the edge ej+i that minimizes 
r on AGj. We take Ei+\ = E i U {ej + i} and let Gj+i be the graph induced 
by the edge set Ei + \. The graph Gj is called the invaded region at time i, 
and the graph S = U°^ Gj is called the invasion percolation cluster (IPC). 
Let Eoo = Ui^o-^- 
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Since we would like to compare Bernoulli percolation to the invasion, we 
use a well-known analogous definition of Bernoulli percolation that makes 
the coupling of the two models immediate. For any p G [0, 1], we say that 
an edge e G E rf is p-open if r e < p. It is obvious that the resulting random 
graph of p-open edges has the same distribution as the one obtained by 
declaring each edge of E d open with probability p and closed with proba- 
bility 1 — p, independently of the states of all other edges. The percolation 
probability 9{p) is the probability that the origin is in the infinite cluster 
of p-open edges. There is a critical probability p c = inf{p : 9{p) > 0} G (0, 1). 
For general background on Bernoulli percolation, we refer the reader to [5]. 

It was shown in [4] that for all p> p c , the invasion intersects the infinite 
p-open cluster with probability 1. In the case d = 2, this result immediately 
follows from the Russo-Seymour- Welsh theorem (see Section 11.7 in [5]). 
Furthermore, the definition of the invasion mechanism implies that if the 
invasion reaches the p-open infinite cluster for some p, then it will never 
leave this cluster. Combining these facts yields that if is the edge added 
at time i, then limsupj^oo r e; = p c . From now on, we consider only d = 2. In 
this case, it is well known that 9{p c ) = 0, which implies that for every t > 0, 
there is an edge e(t) such that e(t) is invaded after step t and T e M > p c . 
The last two results give that f\ = max{r e : e G £"00} exists and is greater 
than p c . Let e\ denote the edge at which the maximum value of r is taken 
and assume that e\ is invaded at step i\ + 1. Following the terminology 
of [15], we call the graph the first pond of the invasion and denote it 
V\. The edge e\ is called the first outlet. The second pond of the invasion 
is defined similarly. Note that the same argument as above implies that 
T2 = max{r ei : a G -Eoo>£ > h} exists and is greater than p c . If we assume 
that T2 is taken on the edge e.2 at step 12 + 1 , we call the graph G7j 2 \ G^ the 
second pond of the invasion and denote it V2. The further ponds Vk can be 
defined analogously. 

The following interpretation gives a natural meaning to the ponds. Con- 
sider an infinite piece of land divided into square parcels. These parcels are 
separated by dikes whose heights are given by the values of independent 
random variables, uniformly distributed on [0, 1]. One of the parcels, called 
the parcel of the origin, contains an infinite source of water. First, the wa- 
ter level in the parcel of the origin rises until it reaches the height of the 
lowest adjacent dike and then spills over into the parcel on the other side 
of this dike. Next, the water level rises in both parcels until it reaches the 
height of the lowest dike on the boundary of the union of the two parcels, at 
which time a new parcel floods. The process continues indefinitely and, as 
time approaches infinity, an infinite region of land will flood. Consider the 
dual lattice of Z 2 , each dual edge having the r value of its corresponding 
edge in the original lattice, identifying the dual edges with the dikes and the 
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origin with the source of water. Each vertex of Z 2 corresponds to exactly 
one parcel of land. It is evident from the invasion mechanism and from the 
way the flood spreads on the land that a parcel is flooded if and only if the 
corresponding vertex of 1? is invaded. We now explain the meaning of the 
first pond in the flood setting. At step i±, when the first outlet is invaded, 
the minimal r value on the boundary of Gi 1 is that of e\. However, this is 
the edge with the largest r value ever added to the invasion. This means 
that the invasion will never return to G^, that is, no edge on AG^, other 
than ei, will be invaded. Therefore, after some time, all water will flow over 
the dike corresponding to e\ and the water level in each parcel of the first 
pond will be constant and equal to t\. The same argument shows that after 
some time, the water level in the second pond will become, and remain, T2, 
and so on. 

Now that our model is defined, we review a few results that established 
connections between the invasion and the critical percolation models. To 
the best of our knowledge, the first paper with mathematically rigorous 
results in this area was [4], where it was shown, among other things, that 
the empirical distribution of the r value of the invaded edges converges to the 
uniform distribution on [0,p c ]. Results on the fractal nature of the invaded 
region were also obtained in [4]. The authors showed that the region has 
zero volume fraction, given that there is no percolation at criticality, and 
that it has boundary-to- volume ratio (1 — p c )/Pc- This corresponds to the 
asymptotic boundary-to-volume ratio for large critical clusters (see [10] and 
[14]). The above results indicate that a large proportion of the edges in the 
IPC belong to big p c -open clusters. 

An object that turns out to be closely related to the invaded region is 
the incipient infinite cluster (IIC). Loosely speaking, one can say that the 
IIC is the "infinite open cluster at criticality." The IIC can be constructed 
by conditioning on the origin being connected to a site at distance n from 
the origin in critical percolation and by considering the cluster of the origin. 
If we let n — > oo, an infinite cluster is obtained and this cluster is called 
the incipient infinite cluster. (Later in this paper, we will give the precise 
definition. For detailed results on the IIC, we refer the reader to [8].) Let S n 
be the number of invaded sites within a distance of at most n from the origin. 
The scaling of the moments of S n as n goes to infinity was obtained in [6] 
and [22] , and it turned out to coincide with the scaling of the corresponding 
moments for the IIC. Another similarity established in [6] is concerned with 
the invasion picture far away from the origin: the invasion measure was 
shown to be locally the same as the IIC measure. 

The diameter and volume of the first pond of the invasion were studied 
in [18, 19]. It was shown that the decay rates of their distributions coincide, 
respectively, with the decay rates of the distributions of the diameter and 
the volume of the critical cluster of the origin in Bernoulli percolation. 
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To the best of knowledge, the only paper to date concerned with the 
differences between the invasion model and critical percolation is [1]. The 
authors consider invasion percolation on regular trees. The scaling behavior 
of the r-point function and the volume of the invaded region at and below a 
given height can be explicitly computed. It is found that while the power laws 
of the scaling are the same for the invaded region and for the incipient infinite 
cluster, the scaling functions differ and, consequently, the two clusters behave 
differently. In fact, their laws are found to be mutually singular. Even though 
the arguments of [1] do not work for invasion in the plane, their results give 
a strong indication that, in spite of the presence of many similarities, the 
two objects are indeed different. 

In this paper, we compare connectivity properties of the origin's invaded 
region to those of the critical percolation cluster of the origin and the IIC. 
In Theorems 1.1 and 1.2, we give the asymptotic behavior for the /c-point 
function of the first pond. We continue to study the relation between the IPC 
and large p c -open clusters in Theorems 1.3 and 1.4. We show that, for any 
K and N, there are infinitely many ponds that contain at least K disjoint 
p c -open clusters of size at least N . We also show that, provided the radius 
of the first pond is larger than TV, the first pond contains at least K disjoint 
p c -open clusters of size at least N with probability bounded from below by a 
positive constant independent of N . For k > 1, we compute the exact decay 
rate of the distribution of the radius of the kth pond in Theorem 1.5. Unlike 
the decay rate of the distribution of the radius of the first pond [18], it is 
strictly different from that of the radius of the critical cluster of the origin. 
Finally, in Theorem 1.8, we show that the IPC measure and the IIC measure 
are mutually singular. 

1.2. Notation. In this section, we set out most of the notation and defi- 
nitions used in the paper. 

For a G R, we write \a\ for the absolute value of a and, for a site x = 
(x\,X2) G Z 2 , we write |x| for max(|xi|, |x2|). For n > and x G Z 2 , let 
B(x, n) = {y G 1? : \y — x\ < n} and dB(x, n) = {y G 1? : \y — x\ = n}. We 
write B(n) for B(0, n) and dB(n) for dB(0,n). For m <n and x G Z 2 , we 
define the annulus Ann(x; m, n) = B(x, n) \ B(x, m). We write Ann(m, n) for 
Ann(0; m, n). 

We consider the square lattice (Z 2 , E 2 ) , where E 2 = {(x,y) GZ 2 xZ 2 :|x — 
y\ = l}. Let (Z 2 )* = (l/2,l/2)+Z 2 and (E 2 )* = (1/2, 1/2) + E 2 be the ver- 
tices and the edges of the dual lattice. For xeZ 2 , we write x* for x + 
(1/2, 1/2). For an edge e G E 2 , we denote its ends, left (resp., right) or bottom 
(resp., top), by e x , e y G 1? . The edge e* = (e x + (1/2, 1/2), e y - (1/2, 1/2)) is 
called the dual edge to e. Its ends, bottom (resp., top) or left (resp., right), 
are denoted by e* and e* Note that, in general, e* and e* are not the same 
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as {e x )* and (e y )* . For a subset K C Z 2 , let K* = (1/2, 1/2) +K. We say that 
an edge e € E 2 is in /C C Z 2 if both of its ends are in /C. 

Let (r e ) ee E2 be independent random variables, uniformly distributed on 
[0, 1], indexed by edges. We call r e the weight of an edge e. We define the 
weight of an edge e* as r e * = r e . We denote the underlying probability mea- 
sure by P and the space of configurations by ([0, 1] E where T is a 
natural c-field on [0, 1] E . We say that an edge e is p-open if r e < p and 
p- closed if r e > p. An edge e* is p-open if e is p-open and it is p- closed if 
e is p-closed. The event that two sets of sites K,\,K,2 C Z 2 are connected 
by a p-open path is denoted by JC\ < — > K-2 and the event that two sets of 
sites Kl,K,2 C (Z 2 )* are connected by a p-closed path in the dual lattice is 

denoted by K\ > . 

For positive integers m<n, k and p € [0,1], let A niP be the event that 
there is a p-open circuit around the origin of diameter at least n and let 
B n ,p be the event that there is a p-closed circuit around the origin in the 
dual lattice of diameter at least n. Let A m ^ n ^ p be the event that there is a 
p-open circuit around the origin in the annulus Ann(m, n) and let B m ^ n ^ v be 
the event that there is a p-closed circuit around the origin in the annulus 
Ann(m,n)*. Let A^ nnp be the event that there are k disjoint p-open paths 
connecting B(m) to dB(n). 

For p € [0,1], we consider a probability space (fip^pjPp), where U p = 
{0, 1} E , J- p is the cr-field generated by the finite-dimensional cylinders of 
Q p and P p is a product measure on (f2 p ,.Fp), P p = IleeE 2 Me, where fj, e is 
given by /i e (w e = 1) = 1 — p e (w e = 0) = p for vectors (w e ) eg]E 2 G fi p . We say 
that an edge e is open or occupied if cj e = 1, and e is closed or vacant if 
a; e = 0. We say that an edge e* is open or occupied if e is open, and it is 
closed or vacant if e is closed. The event that two sets of sites K,\,K,2 C Z 2 
are connected by an open path is denoted by K,\ <-> /C2 and the event that 
two sets of sites /CJ , /C2 C Z 2 are connected by a closed path in the dual 
lattice is denoted by JC \ A /Q . 

For positive integers m < n and k, let A n be the event that there is an 
occupied circuit around the origin of diameter at least n and let B n be the 
event that there is a vacant circuit around the origin in the dual lattice of 
diameter at least n. Let A m ^ n be the event that there is an occupied circuit 
around the origin in the annulus Ann(m,n) and let B m ^ n be the event that 
there is an vacant circuit around the origin in the annulus Ann(m,n)*. Let 
n be the event that there are k disjoint occupied paths connecting B(m) 
to dB(n). 

For two functions g and h from a set X to R, we write g(z) >c h(z) to 
indicate that g{z)/h{z) is bounded away from and 00, uniformly in z € X . 
Throughout this paper, we write "log" for log 2 . We also write P cr for P Pc . 
All of the constants (Cj) in the proofs are strictly positive and finite. Their 
exact values may be different from proof to proof. 
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1.3. Main results. 

1.3.1. Probability for k points in the first pond. 

Theorem 1.1. Let C(0) be the cluster of the origin in Bernoulli bond 
percolation. For any k > 0, 

(l) F(x 1 ,...,x k eV l )~F cr (x 1 ,...,x k eC{0)), Xl ,...,x k eZ 2 . 

Remark 1. The lower bound follows from the observation that the p c - 
open cluster of the origin is a subset of V\. 

The reader may ask whether there is a universal constant c such that, for 
all k > 1 and x\,...,x k GZ 2 , 

P(xi, . . . ,x k E t>i) < cP cr (xi, . . . ,x k e C(0)). 

In the next theorem, we show that the answer to the above question is 
negative. 

Theorem 1.2. 

P(B(n) C Vi) 

im = no 

n-oo F cr (B(n) C C(0)) 

1.3.2. Ponds and p c -open clusters. We now state two theorems which 
say that invasion ponds can contain several large p c -open clusters. Let K > 
2, N > 1, and let U(m,K,N) be the event that the mth pond contains at 
least if disjoint p c -open clusters of size at least N. 

Theorem 1.3. With probability one, there exist infinitely many values 
of m for which U(m, K, N) holds. 

Theorem 1.4. There exists e > 0, independent of N but dependent on 
K, such that 

WQA{l,K,N)\R x >N)>e, 
where R\ is the radius of the first pond. 

1.3.3. Radii of the ponds. We define Rj to be the radius of the graph 
Gi } , that is, Rj = max{|x| :x € Gj-}. We refer the reader to Section 1.1 for 
the definitions of ij and Gi j . In the next theorem, we give the asymptotics 
for the radii Rj. 
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Theorem 1.5. For any k>l, 
(2) F(R k > n) x (logri^P^O <-► dB{n)). 

Remark 2. Let {0 ^ k dB(n)} be the event that there is a path con- 
necting the origin to the boundary of B(n) such that at most k of its edges 
are closed. If this event holds, then we say that the origin is connected to 
dB(n) by an open path with k defects. It is a consequence of the Russo- 
Seymour- Welsh (RSW) theorem (see [17], Proposition 18) that 

P cr (0 ^ k dB{n)) x (logra) fc P cr (0 <-> dB(n)). 

Therefore, Theorem 1.5 implies that, for any k > 1, 

P(fl fc >n)xP cr (0^fc_i5J3(n)). 

Remark 3. For k = 1, the statement (2) follows from Theorem 1 in 
[18]. Note that in the case k = 1, the lower bound immediately follows from 
the fact that C(0) C V\, where C(0) is the p c -open cluster of the origin for 
Bernoulli bond percolation. However, in the case k>2, the lower bound is 
not trivial. 

Let R k be the diameter of the kth pond, R k = max{|x — y\:x,y £ V k }. 
Note that (Rk) are related to (Rk) via the simple inequalities i?i < R\ < 2R\ 
and Rk — Rk-i — 1 < Rk < 2Rk for k > 2. The next theorem immediately 
follows from Theorem 1.5 and the fact that W cr (0 <-> dB(n)) x P C r(0 <-> 
dB(2n)). 

Theorem 1.6. For every k>l, 

F(R k > n) x (logn) fe ~ 1 P cr (0 ^ dB(n)). 

1.3.4. Mutual singularity of IPC and IIC. First, we recall the definition 
of the incipient infinite cluster from [8]. It is shown in [8] that the limit 

u(E)= lim P cr (E\0^dB(N)) 

N~*oo 

exists for any event E that depends on the state of finitely many edges in 
E 2 . The unique extension of v to a probability measure on configurations of 
open and closed edges exists. Under this measure, the open cluster of the 
origin is a.s. infinite. It is called the incipient infinite cluster (IIC). Recall 
the definition of the IPC S from Section 1.1. The next statement is [6], 
Theorem 3. 
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Theorem 1.7. For any finite K. C E 2 and x G Z 2 , Ze£ /C(x) = x + /C C E 2 , 
^ = {/Cc S} and E' K = {K,C (7(0)}. T/ien, 

lim P(£je (x) |a:€ 5) =!/(£&). 

The above theorem says that, asymptotically, the distribution of invaded 
edges near x is given by the IIC measure. In this paper, we show that, 
globally, the IPC measure and the IIC measure are entirely different. 

Theorem 1.8. The laws of IPC and IIC are mutually singular. 

1.4. Structure of the paper. We define the correlation length and state 
some of its properties in Section 2. We prove Theorem 1.1 in Section 3 and 
Theorem 1.2 in Section 4. The proofs of Theorems 1.3 and 1.4 are given in 
Section 5. In Section 6, we prove Theorem 1.5. Theorem 1.8 is proved in 
Section 7. After Sections 1 and 2, the remainder of the paper may be read 
in any order. For the notation in Sections 3-7, we refer the reader to Section 
1.2. 

2. Correlation length and preliminary results. In this section, we de- 
fine the correlation length that will play a crucial role in our proofs. The 
correlation length was introduced in [3] and further studied in [9]. 

2.1. Correlation length. For positive integers m,n and p £ (p c , 1], let 
a(n,m,p) =P p (there is an open horizontal crossing of [0,n] x [0,m]). 
Given e > 0, we define 

(3) L(p, e) = min{n : a(n, n,p) > 1 — e}. 

L(p,e) is called the finite-size scaling correlation length and it is known that 
L(p,e) scales like the usual correlation length (see [9]). It was also shown 
in [9] that the scaling of L(p, e) is independent of e, given that it is small 
enough, that is, there exists £q > such that for all <£i,£2 <£o 5 we have 
L(p, £i) x L(p, £2). For simplicity, we will write L(p) = L(p, £q) for the entire 
paper. We also define 

p n = sup{p:L(p) > n}. 

It is easy to see that L(p) — > 00 as p — > p c and L(p) = 1 for p close to 1. In 
particular, the probability p n is well defined. It is clear from the definitions 
of L(p) and p n , and from the RSW theorem, that for positive integers k 
and I, there exists 5k,i > such that for any positive integer n and for all 

PG [PcPn], 

P p (there is an open horizontal crossing of [0,/cri] x [0, In)) > 5kj 
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and 

P p (there is a closed horizontal dual crossing of ([0, kn] x [0, in))*) > 5^ \. 

By the FKG inequality and a standard gluing argument [5], Section 11.7, 
we get that, for positive integers n and k>2, and for all p € [PoPn], 

P p (Ann(n, kn) contains an open circuit around the origin) > (<52fc,fc-i) 4 

and 

P p (Ann(n, kn)* contains a closed dual circuit around the origin) > (52k,k-l) • 

2.2. Preliminary results. For any positive I, we define log^ ^ / = I and 
log^ I = log(log^' -1 ^ I) for all j > 1, provided the right-hand side is well 
defined. For I > 10, let 

(4) log* I = min{j > : log^ / is well defined and log^ I < 10}. 

Our choice of the constant 10 is quite arbitrary; we could take any other 
large enough positive number instead of 10. For I > 10, let 



(5) Pi{j) 



inf \p > p c : L{p) < 1 X, if j G (0, log* I), 



Pc if J > log* Z, 

1, if ./ = 0. 

The value of C* will be chosen later. Note that there exists a universal 
constant Lo(C*) > 10 such that pi(j) are well defined if I > Lq(C*) and 
nonincreasing in I. The last observation follows from the monotonicity of 
L(p) and the fact that the functions l/log^l are nondecreasing in I for 
j £ (0,log*0 and I > 3. 

We give the following results without proofs: 

1. (Reference [6], (2.10).) There exists a universal constant D\ such that, 
for every I > Lo(C*) and j 6 (0,log* I), 

(6) C* logW I < - \ < D X C* log^ I. 

L (pl(j)) 

2. (Reference [9], Theorem 2.) There exists a constant D2 such that, for all 

P>Pc, 

(7) 0(p) <F P [0~ dB(L{p))] < D 2 F cr [0 <-> dB(L(p))], 

where 6{p) = ¥ p (0 <-> 00) is the percolation function for Bernoulli perco- 
lation. 

3. (Reference [16], Section 4.) There exists a constant D3 such that, for all 
n > 1, 

(8) P Pn (B(n)^oo)>D 3 . 
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4. (Reference [9], (3.61).) There exists a constant D4 such that, for all pos- 
itive integers r < s, 

/q\ F cr (0~dB(s)) fr 

iUj P cr (0~ dB(r)) -^ 4 Vs" 

5. Recall that B n is the event that there is a closed circuit around the origin 
in the dual lattice with diameter at least n. There exist positive constants 
D5 and Dq such that, for all p> p c , 

(10) F p (B n ) < D 6exp {-i} 6 " 

This follows from, for example, [6], (2.6) and (2.8) (see also [17], Lemma 
39 and Remark 40). 

6. (Reference [17], Proposition 34.) Fix e = ((0, 0), (1, 0)) and let A% 2 be the 
event that e x and e y are connected to dB(n) by open paths, and e* and 
e* are connected to dB(n)* by closed paths. Note that these four paths 
are disjoint and alternate. Then, 

(11) (Pn ~ P c)n 2 F cr (A% 2 ) x 1, n > 1. 

3. Proof of Theorem 1.1. Before we prove Theorem 1.1, we give two 
lemmas that will be used in the proof. To simplify the notation, we write 
= xq. For positive integers m <n and i£Z 2 , we define the event 

(12) A mtTl (x) = {there is an open circuit in the annulus Ann(x; m,n)}. 

Lemma 3.1. Given a set of vertices {x\, . . . , x^} G Z 2 , let mi = min{|xj — 
Xj\-0<j<k,j^i}, where xq = and let m = minjmj : < i < k}. Further- 
more, assume m = m k - There then exists a constant C\, independent of k, 
such that for all p> p c , the probability 

P p (xi «-> 00, . . . , Xk *-* 00, <-> 00) 

is bounded from above by 

Ci~P p (xk <-> dB(xf c ,m))F p (xi <-> 00, . . . , x^-i <-> 00, <-> 00). 

Proof. The statement is trivial if m < 4, so we assume that m > 4. By 
the RSW theorem, there is a constant C 2 independent of k and m such that, 
for &\\p>p c , Pp(^[ m / 4 ],[ m /2](^fc)) > 1/C 2 and hence 1 < C 2 Pp(^[ m / 4 ],[m/2] (^fc))- 
The FKG inequality gives 

P p (xi «-> 00, . . . , <-> 00, <-> 00) 

(13) < C7 2 Pp(A[ m /4] i [ m / 2 ](xfc))Pp(xi <-> 00, . . . ,x fc <-> oo,0^ 00) 

< C2Pp(A[ m /4] i [. m / 2 ](xfc),Xi <-> OO, . . . , X k <-» OO,0 <->• OO). 

The event on the right-hand side of (13) implies the following two events: 
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1. {x k ^dB(x k , [m/4])}; 

2. {xi <-> oo, . . . , <-> oo, <-> oo outside B(x k , [m/4])}. 

These two events are independent and therefore the right-hand side of (13) 
is bounded from above by 

C 2 F p (x k ^dB(x k ,[m/4])) 

x ¥ p (xi <-> oo, . . . , <-> oo, <-> oo outside B{x k , [m/4])) 

< C2P p (xfc <-> dB(x k , [m/4\))¥ p (xi oo, . . . , Xfe-i oo, <-> oo), 

where the last inequality follows from monotonicity. Finally, it follows from 
the FKG inequality, RSW theorem and a standard gluing argument [5], 
Section 11.7, that ¥ p (x k <-> dB(x k , [m/4])) x P p (x k <-> dB(x k , m)) uniformly 
in p > p c . □ 

We recall the definition of the probabilities (p n (j)) in (5). We also re- 
call that these probabilities are well defined if n > Lo(C*), where C* is the 
constant from (5). Later, we choose C* to be sufficiently large. 

Lemma 3.2. Given a set of vertices {x\, . . . ,x k } € Z 2 , let n = max{|xj — 
Xj\:i,j = 0,...,k}, where xq = 0. Furthermore, assume that n > Lo(C*). 
There is then a universal constant C3 such that, for all j € (0,log*n), 

^p n (j)( x l ^ 00,. . . ,x fe ^ oo,0^ 00) 

< (C 3 logW n) (fc+1)/2 P cr (x!, . . . , x fc € C(0)). 

Proof. We will use induction in A:. First, we consider the case k = 1. To 
simplify our notation, we write £1 = x. Note that, now, \x\ = n = m, where 
m is defined as in Lemma 3.1. From Lemma 3.1, it follows that 

(15) f Pn (j)(x oo,0 <-» 00) < Ci0(p n (j))P Pn(i) (O <- dB(n)). 
Since L(p n (j)) < n, we obtain 

P Pn(j) (0 <-> 0B(n)) < F Mj) (0 <-> dB(L( Pn (j)))). 
Combined with (6), (7) and (9), the above inequality gives 

C x 0(jp n {j))F pn{j) (Q ~ dB(n)) < C 4 P cr (0 ~ dB(L( Pn (j)))) 2 

< C 5 — " P cr (0 ~ dB{n)f < C 6 logW nP cr (0 <-► dB{n)f. 

L (Pn{j)) 

The RSW theorem and the gluing argument show (see, e.g., [7], (4)) that 

(16) P cr (0 «-» dB{n)f < C 7 P cr {x G (7(0)) 



INVASION PERCOLATION IN 2D 



13 



for some constant C7. In particular, (14) follows for k = 1. 

The general case is more involved. We assume that Lemma 3.2 is proved 
for any set of vertices {y\, . . . , y k ~i} G 1? . Then, for a set of vertices {x±, . . . , x k } G 
7?, we define m as in Lemma 3.1 and assume that m = m k = min{|xj — 
x k \ : i < k}. We also define n\ = max{|xj — xj \ : i, j = 0, . . . , k — 1}, with xq = 
0. Then, by the induction hypothesis, 

^Pni (j)( X l °0,.. -,Xk-l <-> OO,0 <-> OO) 

(17) 

< (C7 3 logW m^P^i, . . . , G C(0)). 

Since n x < n and m < n, we get p n (j) < p m (j) and p n (j) < p ni (j) (see Sec- 
tion 2). Therefore, 

^Wo')^ 1 ^ 00 ' • • • > Xk ^ ^ 00 ) 

< CiPp m(j -)(a; fe <-> 95(x fc ,m))P pni(j )(a;i <->oo,.. .,x k -i 00, <-> 00) 

< CiP^^^fc 9J3(x fc ,m))(C 3 log^ nx) fc/2 P cr (xi, . . . ,x fc _i G C(0)) 

< (C7 8 log^m) 1/2 P cr (x fc ^5S(x fe ,m)) 

x (C 3 log® n 1 ) fe/2 P CT (x 1 , . . . , x fe _x G C(0)) 

< Cj /2 C* /2 (log^ n) (fc+1)/2 P cr (x fc <- 3B(x k , m))P cr (x!, . . . , x fc _x G C(0)), 

where the first inequality follows from Lemma 3.1 and monotonicity, the 
second inequality follows from (17) and the third inequality follows from (6) 
and (9). Note that Cs is independent of k. It now suffices to show that there 
is a universal constant C9 such that 

P cr (x k <-> dB(x k ,m))F cr (xi, . . .,x k -i G C(0)) 

(18) 

<C 9 F cr ( Xl ,...,x k eC(0)). 

Assume that (18) is proved. We can then take C3 = max-fC^Cy, C%Cg}. The 
argument above shows that we can proceed to the next k using this value of 
C3. We now show (18). We take X{ such that m = \x k — Note that this 
vertex may be the origin. We know that at least one such vertex exists. Recall 
the definition of events A m ^ n (x) from (12). By the RSW theorem, there is 
a constant C w such that 1 < C w F cr (A[ m / 2 ] :m ( x i)', A [rn/2],m( x k))- Using the 
FKG inequality, we get 

F cr (x k <-> dB(x k ,m))F cr (xi, . . .,x k -i G C(0)) 

< C w F cr (A [m /2i m {xi);A [m/2 ], m (xk))^cr{x k <-> dB(x k ,mj) 
x P CP (xi,...,x fc _i EC(0)) 

< CioP C r(^4[m/2],m( a:; i); ^4 

[m/2],m 

x k ^dB(x k ,m);x 1 ,...,x k -i G C(0)). 
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We show that the event 

{A[ m /2],m( x iy> A [m/2],m( x k)]Xk <~> 9B(x k , m); X\, . . . € C(0)} 

implies the event {xj <->• x k ;x\, . . . , x k -i € C(0)}. Indeed, it follows from sim- 
ple observations: 

1. Since the events {xk <-> dB(xf t ,m)} and A\ m j2\ m {x] t ) hold, x k is con- 
nected to the circuit lying in the annulus Ann(x£,; [m/2],m). 

2. Since the distance between Xj and xj, is m, the boxes B(xi, [m/2] + 1) 
and B(xk, [m/2] + 1) intersect. This implies that the circuits in the annuli 
Ann (x j.; [m/2], m) and Ann(xj; [m/2], m) intersect. 

3. Recall that m is the minimal distance in the graph with vertex set 
{0, x\, . . . , Xk}- Since > 2 and {xi, . . . , x k -x S C(0)}> there is a vertex 
#j 7^ Xk (it may be the origin) such that Xj ^ B(xi,m — 1) and x_j is con- 
nected to Xj. The last observation implies that x« is connected to the 
circuit lying in Ann(xj; [m/2],m) and hence also to x k . 

This proves (18). □ 

Proof of Theorem 1.1. For {x l5 . . . ,x/J G Z 2 , we define, as in Lemma 
3.2, n = max{|xj — Xj \ : i, j = 0, . . . , k}. If n < Lq(C*)> then P cr (xi, . . . , Xk £ 
C(0)) > const (C*). Theorem 1.1 immediately follows since P(xi,...,Xfe € 
Vi) < 1. We can therefore assume that n > Lo(C*). In particular, the proba- 
bilities p n (j) are well defined. The rest of the proof is similar to the proof of 
Theorem 1 in [18]. Recall that f\ is the value of the outlet of the first pond. 
We decompose the event {xi, . . .x k € V±} according to the value of f\. We 
write 

log* n 

(19) F( Xl ,...,x k eVi)= £ ¥{x 1 ,...,x k £Vi,f 1 £[p n {j), Pn {j -1))). 
Note that, for any p> p c , 

(a) if fi < p, then any invaded site is in the infinite p-open cluster; 

(b) if a given set of vertices {xi, . . .,x k } is in the first pond, n is defined 
as in Lemma 3.2 and f\ > p, then there is a p-closed circuit around the 
origin with diameter at least n. 

We recall the definition of the event 

B n ,p = {3p-closed circuit around in the dual with diameter at least n}. 

We conclude that the probability P(xi, . . . ,x k 6 V\,t\ 6 [p n (j),Pn(j — 1))) is 
bounded from above by 

(20) P(xi < — > oo,...,x k < — > oo,0 < — > oo\B n ^ n{j) ). 
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The FKG inequality implies that the probability (20) is not bigger than 

F Pn (j-i)( x i °°> • • • > x k <-> oo, oo)P( J B n)Pn(i) ) 
(21) _ 

< Cn(log (i ^n) 12 Pp n (j-i)(a;i ^oo,...,x fc ^oo,0^oo), 

where we use (6) and (10) to bound the probability of B n>p u\ by 

Cii(log (j ~ 1) n)- Cl2 . The constant C 12 can be made arbitrarily large pro- 
vided that C* is made large enough. We consider bounds for (21) separately 
for j = 1 and for j > 1. If j > 1, we use Lemma 3.2 to bound (21) by 

daOog^ 1 ) nf +1 V 2 - C ™F cr (xi, ■ ■ ■ ,x k G C(0)). 
If j = 1, we bound (21) by 
Cnn- Cl2 < Ci 4 n~ 1/2 P cr (0 <-> dB(n)) 2k < C 15 n~ 1/2 F cr (x u ...,x k € (7(0)). 

The first inequality holds for Ci 2 > ft + 1/2 since P cr (0 <-> dB{n)) > in" 1 / 2 
(see [5], (11.90)). The last inequality follows from (16), applied k times, 
and the FKG inequality. Therefore, for all j, if C\i > k + 1/2, then (21) is 
bounded by 

CiaQogG- 1 ) n)- 1/2 P cr (xx, ...,x k € C(0)). 
We plug this bound into (19): 

log* n 



(xi, . .. ,ijfc G Vi) < Ci 6 P cr (xi, . .. G C(0)) ^ (!og G 



^0-1)^-1/2 

<C 17 P cr (x 1 ,...,x fc GC(0)). 
The last inequality follows from the fact that 

log* n 

sup (log^ -1 -* n) -1 / 2 < oo 
n>10 ^ 

(see, e.g., [6], (2.26)). □ 

4. Proof of Theorem 1.2. In this section, we prove that 



hm -v-» cy i) =qo 
n-oo F cr (B(n) C C(0)) 

By RSW arguments [5], Section 11.7, the denominator is at most equal to 

CiFcr(B(n) C C(0) in B{2n)) 
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for some C\ > 0. Recall that p n = sup{p: L(p) > n}. We can bound the 
numerator from below: it is at least equal to 

F Pn( B ( n ) C C(0) in B(2n) n 3 closed circuit around B(2n)) 

= F Pn (B(n) C (7(0) in B(2n))F Pn (3 closed circuit around B(2n)). 

By the definition of L(p), there exists C*2 > such that this probability is at 
least 

C 2 ¥ Pn (B(n)cC(0) in B(2n)). 

Therefore, to prove Theorem 1.2, it suffices to show that 

(99) v F Pn (B(n)cC(0) in B{2n)) 

1 j ™ P cr (B(ri) C C(0) in B(2n)) °°' 

For this, we use Russo's formula [5] (the definition of pivotal edges is also 
given in [5]). Let T n be the event which appears both in the numerator and 
in the denominator of (22). Let p € [e, 1 — e] for some e < \ and, for any 
vertex v, let E v be the set of edges incident to v. We see that 

|-P p (r n ) = £> p (e is pivotal for T n ) 

> — X] p p( e is P ivotal for r n ;r n ) 

P veB(n) e£E v 

> — X! ^(^e G -El. pivotal for T n ; T n ) 

P «eB(n) 

>J_ £ min(p,l-p) 4 P p (r n ) 

^ v£B(ri) 



In particular, 



> c 3 n 2 p p (r n ). 

P P „(r n )>p cr (r n )e C4 " 2 ^-^) 



for some C4 > 0. It follows from (11) and the fact that 9(p c ) = that n 2 (p n — 
p c ) — » 00. This completes the proof. 

5. Proofs of Theorems 1.3 and 1.4. First, we prove two lemmas (see 
Section 1.2 for the definitions). 

Lemma 5.1. For each k>2, there exists such that, for all n, 
where — > as k — > 00 . 
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Proof. Recall that B n ^ n = {there is a closed circuit in Ann(n, 2n)*}. 
Pick c > such that, for all N > 1, 

Wcr(B N ,2N) > c. 

We split the annulus Ann(n, kn) into [log A;] disjoint annuli Ann(2 l n, 2 l+1 n): 

P« fcniPc )<(l-c) 1 °s fc - 1 . 
This completes the proof. □ 

Lemma 5.2. There exists C\ > such that for all N and k, 

Proof. By RSW arguments, there exists C 2 > such that for all N and 

k, 

P(Av,27V,p c n AkN,2kN,p c ) > C%- 

If follows from (8) that there exists C3 > such that for all N and k, 

The FKG inequality gives the result. □ 
We now prove the theorems. 

Proof of Theorem 1.3. We prove the theorem for K = 2. For other 
values of K, the proof is similar. Let D(k,N) = A^ t 2N,p c H Aw,2fc7V,p c ^ 
A\f 2 kNp N an d pick Ci from Lemma 5.2. Fix such that the constant c&/ 2 
from Lemma 5.1 satisfies Cfc/ 2 < 2 • ^ f° nows that 

pp(fc,iV)n{A^ fcJViPc } c )>^. 

For any k>2, there exists C4 = C±{k) such that for all iV, 

P(-B2fcAr,4fc7V, PJV ) > C*4- 

Therefore, by independence, 

P(£>(*. n {Al NjkN:Pc } c n i? 2fc 7v,4 fc iv, PJ v) > ^ > 0. 

This statement, along with the Borel-Cantelli lemma, gives the theorem. 

□ 

Proof of Theorem 1.4. Let A\ p = {0 <-> dB{n) by a p-open path}. 
We first note that [18] gives a constant C5 > such that for all TV, 

P(i?i>iV)<C5P(A 27ViPc ). 
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It is obvious thatP( J Ri > NC\U(1, 2, N)) > T(A\ N ^ pc nW(l, 2, N)). Therefore, 
it suffices to show that there is an e > such that for all N, 

F(U(l,2,N)\Al N>p J>s. 

The rest of the proof is almost the same as the proof of Theorem 1.3. Let 
D{k,N) be as in the proof of Theorem 1.3. Pick C\ from Lemma 5.2. By 
the FKG inequality, we see that 

2N,pJ- 

By independence and Lemma 5.1, we may fix k such that for all N, 

nAi N>Pc nAi N>kN>p j < c k/2 nA N , Vc ) < Y p (^k Pc )- 

For any k>2, there exists C4 = C^(k) such that for all N, 
Independence now gives us 

nAN, Pc n D(k, N) n {Al N ^ Pc y n B 2kNAkN)PN ) > 9^±^ A \ Npc ). 

This concludes the proof. □ 
6. Proof of Theorem 1.5. 



6.1. Upper bound. We give the proof for k = 2. The case k = 1 is consid- 
ered in [18] and the proof for k > 3 is similar to the proof for k = 2. 
We fix n and divide the box B(n) into [logn] + 1 annuli. We write 



[logn]+l 

(23) ¥(R 2 >n)=¥(R 1 >n)+ ^ ¥[R 2 >n,R 1 e 

k=l 



n n 



2 k ' 2 k ~ 



Since [18], Theorem 1, 

P(i?i > n) < CiP cr (0 ^ dB(n)) < CilognP cr (0 <-> dB(n)), 

it remains to bound the typical term of the sum on the right-hand side of 

(23) . It is sufficient to show that there exists a constant C 2 such that, for 
any m £ [0,n/2], 

(24) F(R 2 >n;Ri G [m, 2m]) < C 2 P cr (0 dB{n)). 

We only consider the case m > Lq(C*). The proof for m < Lq(C*) is similar 
to the proof for m > Lo(C*), but much simpler. We omit the details. We 
now assume that m > Lo(C*)- I n particular, the probabilities (p m {i)) and 
(p n (j)) are well defined. 
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We decompose the event on the left-hand side according to the r value 
of the first and the second outlet. The probability P(i?2 > n;R\ £ [m,2m]) 
is bounded from above by 

log* m log* n 

E P (^2 >n;RxE[m, 2m] ; 

1=1 7 = 1 

(25) 

Tl G \PmW,PmV ~ l)];r 2 € [Pn{j),Pn{j ~ !)])• 

Note that if the event {i?i > m;fi 6 [p TO (^);Pm(« — 1)]} occurs, then: 

- there is a p m (i — l)-open path from the origin to infinity; 

- the origin is surrounded by a p m (i)-closed circuit of diameter at least m 
in the dual lattice. 

We also note that if the event {R\ < 2m; R2 > n; T2 G \Pn(j),Pn(j — 1)]} 
occurs, then: 

- there is a p n (j — l)-open path from the box B{2m) to infinity; 

- the origin is surrounded by a p n (j)-closed circuit of diameter at least n in 
the dual lattice. 

From the two observations above, the sum (25) is less than 

log* m log* n 

£ £ P(0 P ^ 1) 5S(m);S(2m) P ^ 1) ^(n);i? m>Pm(i) ;i? niPji0 . ) ). 
i=l j=l 

(26) 

The FKG inequality and the independence of the first two events together 
imply that (26) is not larger than 

log* m log* n 

E E F P^-i)(° " dB(m))P M ^ 1) (B(2m) <- dB{n)) 
i=i 7=1 

(27) 



We use (6) and (10) to bound the probability oi B m>pm ^ by C 3 (log^ _1 ) m)~ Ci , 
where C4 can be made arbitrarily large, provided that C* is made large 
enough. Substitution gives a bound for the last term of (27): 



¥(B mjPmii) ;B n>PriU) ) < min[C 3 (log^-^ m)^ 4 , C 3 (log 



(28) = C 3 max[log (i_1) m, log^'" 1 ' n]"° 4 

< C^log^ 1 ) m)~ C * /2 QogV- 1 '> n)-° 4/2 . 
The RSW theorem and the FKG inequality together imply that 

(29) P p (0 «-» dB(m))F p (B(2m) <-» 0B(n)) < C 5 P p (0 <-> 9J3(n)), 
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uniformly in p >p c - Furthermore, using (6)-(9), we get 

(30) Pfeti-ijlO^SBH) <C 6 (log^rn) 1/2 F cr (0^dB(m)) 
and 

(31) F pn{j _ 1) {B(2m) «-» dB{n)) < C 7 Qog^~^ n) 1/2 F cr (B(2m) <-► dB{n)). 

In the last inequality, we also use (29). We apply the inequalities (28), (29), 
(30) and (31) to (27). We obtain that the probability P(R 2 > n; R\ € [m, 2m]) 
is not larger than 

log* m log* n 

C 8 P cr (0 dB(n)) £ (log (i - 1) m)^ (C4 ~ 1)/2 (log^'- 1 )n)^ (C4 ^ 1)/2 . 

i=l j=l 

We take C* large enough so that C4 is greater than 1. As in (2.26) of [6], it 
is easy to see that there exists a universal constant Cg < 00 such that for all 
n > 10, 

log* n 

6.2. Lower bound. We first give the main idea of the proof. Recall from 
Remark 2 that it is equivalent to prove that P(i?fc > n) > CfcP cr (0 
dB(n)) for some positive constants that do not depend on n. In the 
case k = 1, the event {0 <-^-> <9-B(n)} obviously implies the event {Ri > n}. 
However, for k>2, the event {0 < Pc > k-i dB(n)} does not, in general, imply 
the event {R^ > n}. The weights of some defected edges from the defini- 
tion of the event {0 <^^k-\ dB(n)} can be large enough so that these edges 
are never invaded. We resolve this problem by constructing a subevent of 
the event {0 < c > fc_i dB(n)} which implies the event {R^ > n} and, more- 
over, the probability of this new event is comparable with the probability 
P(0 < Pc > fc_i dB(n)). To construct such an event, we first extend results from 
[9] in Lemmas 6.2 and 6.3 below. We then construct events that will be used 
in the proof of the lower bound in Theorem 1.5 and show that they satisfy 
the desired properties (see, e.g., Corollary 6.2 below). 

We begin with some definitions and lemmas. 

Lemma 6.1 (Generalized FKG). Let £1, . . . ,£„ be i.i.d. real-valued ran- 
dom variables. Let I\,l2,l3 be disjoint subsets o/{l, . . . , n}. Let A\ E o~(£i : i € 
1% U I2) and A<i G <r(£j : i € I2) be increasing in . Let B\ € <r(£j : i G L\ U I3) 
and B2 € <r(£j : i G ^3) 6e decreasing in Then, 



(32) 



p(A 3 nB 3 |^inSi) >p(a 2 )p(s 2 ). 
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Proof. Inequality (32) for P p (rather than P) is given in [9], Lemma 
3, or [17], Lemma 13. The main ingredient of that proof is the Harris-FKG 
inequality for P p (see [5], Theorem 2.4), which is also valid for P (see, e.g., 
[12], Theorem 5.13). Apart from that, the proof of (32) is analogous to the 
proofs of [9], Lemma 3, and [17], Lemma 13, and so we omit it. □ 

Although we will not apply Lemma 6.1 to the following events, they serve 
as simple examples. The events {0 <-^->- dB(n)}, {B(m) dB(n)} are de- 
creasing in (r e ) and the events {0* ^— > dB(n)*}, {B(m)* ^— >dB(n)*} are 
increasing in (r e ). 

Recall that the ends of an edge e € E 2 , left (resp., right) or bottom (resp., 
top), are denoted by e x ,e y £ 1? and the ends of its dual edge e*, bottom 
(resp., top) or left (resp., right), are denoted by e* and e*. We also write 
(1,0) for the edge with ends (0,0), (1,0) G Z 2 . 

Definition 6.1. For any positive integer n, (71,(72 € [0, 1], z € Z 2 and an 
edge e € B(z,n), we define A e (z; n; qi, (72) as the event that there exist four 
disjoint paths P1—P4 such that: 

- P\ and P2 are q\ -open paths in B(z,n) \ {e}, the path Pi connects e x to 
dB(z,n) and the path P2 connects e y to dB(z,n); 

- P3 and P4 are (72-closed paths in B(z,n)* \ {e*}, the path P3 connects e* 
to dB(z,n)* and the path P4 connects e* to dB(z,n)*. 

We write A e (n;qi,q 2 ) for A e (Q;n;qi,q 2 ) and A(n;qi,q 2 ) for A( 1)0 )(n; q 1 , q 2 ). 

For any two positive integers n < N, qi, q 2 6 [0, 1], z£Z 2 , we define 
^4(z; n, A^; (71, (72) as the event that there exist four disjoint paths, two q\- 
open paths in the annulus Ann(z;n,A) from B(z,n) to dB(z,N) and two 
(72-closed paths in the annulus Ann(z; n, A)* from B(z,n)* to dB(z,N)* , 
such that the (71-open paths are separated by the (72-closed paths. We write 
A(n, N; qi,q 2 ) for A(0; n, N; qi,q 2 ). The events A e (n; qi,q 2 ) and A(n, A; gi, g- 2 ) 
are illustrated in Figure 1. 

We will follow the ideas developed in [9]. For that, we need to define 
some subevents of A e (z;n;qi,q2) and A(z; n, A; q±, q 2 ). For n > 1, let U n = 
dB(n) n {x 2 = n}, D n = dB(n) D {x 2 = —n}, R n = dB{n) n {xi = n} and 
L n = dB{n) n {xi = — n} be the sides of the box B(n). Let U n {z) = z + U n , 
D n (z) = z + D n , Rn{z) = z + R n and L n (z) = z + L n be the sides of the box 
P(z,n). 

Definition 6.2. For any positive integer n, q±,q 2 S [0, 1], z£Z 2 and an 
edge e £ B(z,n), we define ^4 e (2;; n; q±, q 2 ) as the event that there exist four 
disjoint paths P\—Pa such that: 
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B(n) 




FlG. 1. Events A e (n;qi,q2) and A(n,N;q 1 ,q2). The solid curves represent q 1 -open paths, 
and the dotted curves represent qi -closed paths. The edge e does not have to be q\ -open or 
qi -closed. 

- P\ and P2 are 91 -open paths in B(z,n) \ {e}, the path Pi connects e x or 
e y to U n {z) and the path P 2 connects the other end of e to D n (z); 

- P3 and P4 are 92-closed paths in B(z,n)* \ {e*}, the path P3 connects e* 
or e* to R n (z)* and the path P4 connects the other end of e* to L n (z)* . 

We define A e (z;n;qi, ■) as the event that there exist two disjoint 91-open 
paths Pi and P2 in B(z,n) \ {e}, the path Pi connects e x or e y to U n (z) 
and the path P2 connects the other end of e to D n (z). 

We write A e (n\ q u q 2 ) for A e (0; n; gi, ? 2 ), ^(™;<?i,92) for A (1)0 )(n; 91,92) 
and we use similar notation for the events A e (z;n;qi, •). 

For any two positive integers n < N, 91,92 £ [0, 1] and z€Z 2 , we define 
j4(z; n, TV; 91, 92) as the event that there exist four disjoint paths P1-P4 such 
that: 

- Pi and P2 are 91-open paths in the annulus Ann(z;n, N), the path Pi 
connects U n (z) to Un(z) and the path P2 connects D n (z) to Dn(z); 

- P3 and P4 are 92-closed paths in the annulus Ann(z; n, N)* , the path P3 
connects R n (z)* to Rm{z)* and the path P4 connects L n (z)* to Ln(z)* . 

We write A(n,N;qi,q2) for ^4(0; n,N; 91, 92). 

We also need to define events similar to the events A in [9], Figure 8. 
For any two positive integers n < N and zGZ 2 , we define U n j^(z) = z + 
[-n, n) x [ti + 1,N], D nyN (z) = z + [-77,77] x [-N, -77 - 1], R U: n(z) = z + 
[77 + 1,N] x [—77, 77] and L n ^(z) = z + [-N, —77 — 1] x [—77,77]. 

Definition 6.3. For any positive integer 77, 91,92 6 [0, 1], z€Z 2 and an 
edge e G B(z, [77/2]), we define A e (z;n; 91, 92) as the event that: 
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B(n/2) 






B(n) 







Fig. 2. Event A e (n;qi,q^). The solid curves represent qi-open paths and the dotted 
curves represent q^-closed paths. The edge e does not have to be qi-open or qi-closed. 

- the event A e (z;n;qi,q 2 ) occurs; 

- the two gi-open paths Pi and P 2 from the definition of A e (z; n; qi,q 2 ) sat- 
isfy Pi nAnn(z; [n/2],n) C U [n/2 ^ n (z) and P 2 n Ann(z; [n/2], n) C D [n/2] ^ n (z); 

- the two ^-closed paths P3 and P4 from the definition of A e (z\ n; qi, q 2 ) 
satisfy P 3 n Ann(z; [n/2],n)* C R^ n / 2 ^ n (z)* and P 4 n Ann(z; [n/2],n)* C 

^[n/2],n( z )*; 

- there exist gi-open horizontal crossings of U\ n i 2 \ n (z) and D\ n / 2 } n (z) and 
there exist g 2 -closed vertical crossings of -^[ n /2],n( z )* an d P[n/2],n( z )*- 

We write A e (n;qi,q 2 ) for A e (0; n; qi, q 2 ) and A(n;gi,g 2 ) for v4 (1)0 )(rj;gi,g 2 ). 

The event A e (n; gi, q 2 ) is illustrated in Figure 2. 

For any positive integers n, iV such that 4n < N, gi,g 2 G [0, 1], z € Z 2 , we 
define A(z;n, N; q\, q 2 ) as the event that: 

- the event A(z;n, N;q\,q 2 ) occurs; 

- the two f/i-open paths Pi and P 2 from the definition of A(z; n, N; q±, q 2 ) 
satisfy P x n Ann(>; n, 2n) C [/„, 2 „(z), Pi n Ann(z; [JV/2], iV) C Efy/g.ivO*), 
P 2 nAnn(2;n,2n)cD n ,2„W and P 2 nAnn(z; [iV/2],^) C D [N/2 ] >N (z); 
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- the two (/2-closed paths P3 and -P4 from the definition of A(z; n,N; 91,92) 
satisfy P 3 nAnn(z;n,2n)* C Rn,2n(z)*, P 3 nAnn(z; [N/2],N)* C R [n/2 ],n(z)* , 
P 4 n Ann(z; n, 2n)* C L ni2n (.z)* and P 4 n Ann(z; [A/2], A)* C £[jv/2],iv0)*; 

- there exist gi-open horizontal crossings of U nt2 n(z), U[n/ 2 ],n(, z )i D n>2n (z) 
and D[ N /2],n( z )> and there exist ^-closed vertical crossings of L nt2n (z)* , 
L{n/2],n( z )* j ^n,2n(^)* and ifyv/^.ivC 2 )*- 

We write A(n, TV; 91, 92) for ^4(0; n, N; q±, 92)- 

Lemma 6.2. For any positive integers n, N such that 4n < N and q±,q2 €E 

(33) P(A(n, A; gi ,q 2 )) x P(A(n, A; gi, g 2 )) ~ P(A(n, A; g 1; g 2 )) 
and 

(34) P(A(AT; q x ,q 2 )) x P(A(A; g 1; g 2 )) x P(I(JV; gi, <&)), 

where the constants in (33) and (34) do not depend on n, N , q± and q 2 - 

Proof. The case gi = g2 is considered in [9], Lemma 4 (see also [17], 
Theorem 11). The proof is based on Lemma 6.1 and the RSW theorem. The 
same proof is valid for general gi and g2- □ 

We need several corollaries of Lemmas 6.1 and 6.2. Their proofs are sim- 
ilar to the proofs for gi = g2 (see, e.g., Corollary 3 and Lemma 6 in [9] or 
Propositions 12 and 17 in [17]). We omit the details. 

Corollary 6.1. 1. For any positive integers a, b and n < N such that 
an < bN, for any q 1 ,q 2 £ \p c ,Pn], 

(35) P(A(n, A; q h q 2 )) X P(A(an, bN; qi ,q 2 )), 
where the constants in (35) only depend on a and b. 

2. For any positive integers n <m < N and qi,q 2 £ \Pc,Pn], 

(36) F(A(n,N;q 1 ,q 2 ))^P(A(n,m;q 1 ,q 2 ))P(A(m,N;q 1 ,q 2 )), 

where the constants in (36) do not depend on n, m, N , q\ and q 2 . 

3. For any positive integer N , qi,q 2 £ \Pc,Pn] and edge e € B([N/2]), 

F(A e (N; gi ,q 2 )) x ¥(A e (N; q 1} q 2 )) X P(i e (A; q u q 2 )) 

(37) 

where the constants in (37) do not depend on N , q\, q 2 and e. 
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The proof of the lower bound in Theorem 1.5 is based on the following 
lemma. 

LEMMA 6.3. For any positive integer N , q±,q2 £ [PcjPn] ande € B{[N/2\), 
(38) F(A e (N; qi ,q 2 ))~F(A(N;p c , Pc )), 
where the constants in (38) do not depend on N , q\, q 2 and e. 

PROOF. The proof for q\ = q 2 is given in [9], Lemma 8, and [17], The- 
orem 27. In this case, the probability measure F can be replaced by the 
probability measure F qi on configurations of open and closed edges. This 
is not the case when q± ^ q 2 , which makes the proof of (38) more involved. 
Note that, by (34) and (37), it is sufficient to show that, for qi,q 2 S \Pc,Pn], 

F(A(N; qi ,q 2 ))~V(A(N;p c , Pc )). 

It is immediate from monotonicity in q\ and q 2 that 

P(A(N;p c ,q 2 )) <F(A(N;q u q 2 )) <F(A(N; qi ,p c )). 

Therefore, it remains to show that there exist constants D\ and D 2 such 
that for all qi,q 2 £ [p c ,Pn], 

¥(A(N; Pc ,q 2 ))>D 1 P(A(N;p c , Pc )) 

and 

F(A(N; qi ,p c ))<D 2 F(A(N;p c ,p c )). 

Since the proofs of the above inequalities are similar, we only prove the first 
inequality. For that, we use a generalization of Russo's formula [5]. We take a 
small 5 > 0. The difference F(A(N;p c ,p)) -F(A(N;p c ,p + 8)) can be written 
as the sum 

5 Yl F(A(N;p c ,-),A e (N;p,-),D e (N;p))+0(5 2 ), 
eeB(N),e^(l,0) 

where D e (N;p) is the event that there exist three p-closed paths P\ — P% 
in B(N)*; the path Pi connects an end of the edge (1,0)* to an end of the 
edge e*; the path P 2 connects the other end of the edge (1,0)* to R* N and 
the path P3 connects the other end of the edge e* to L* N ; or the path P 2 
connects the other end of the edge (1,0)* to L* N and the path P3 connects 
the other end of the edge e* to R* N . Letting 5 tend to 0, we obtain 



(39) — F(A(N;p c ,p)) = - Y,nA(N; Pc ,-),A e (N;p,-),D e (N;p)). 
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We write the right-hand side of (39) as 

[N/2] 

( 40 ) "E E nA(N;Pc,-),MN;p,-),De(N;p)), 

3=1 e: \e x \=j 
N 

(41) E E nMN;Pc,-),A e (N-p,-),D e (N;p)). 

j=[N/2]+le:\e x \=j 

By independence, the sum (40) is bounded from below by 

[N/2] 

- E nM[jm;PcP))nM[^m,N;p c ,p))F(A e (e x ;[j/2];p,p)). 

3=1 e: \e x \=j 

We use (35), the bound (J{e: \e x \ = j} < 16j and the fact that Lemma 6.3 is 
proved for qi = q 2 to bound the above sums from below by 

[N/2] 

-d ]T jF(A(j;p c ,p))F(A(j,N;p c ,p))F(A(j;p c ,p c )) 

3=1 

(42) 

[N/2] 

>-C 2 F(A(N; Pc ,p)) jF(A(j;p c , Pc )), 

3=1 

where the inequality follows from (36). We estimate the sum in (42) using 
the relation 

N 

(43) E^W;^'^)) X N 2 F(A(N; Pc ,p c )). 

3=1 

The relation (43) follows from (36) and the fact that F(A(j,N;p c ,p c )) > 
C?,{j /N) 2 ~ Ci for some positive C3 and C4 that do not depend on j and 
N . This fact follows, for example, from [17], Theorem 24, where the 5-arms 
exponent is computed for site percolation on the triangular lattice. The same 
proof applies to bond percolation on the square lattice. 

Similarly to the proof of [20], Lemma 6.2, the sum (41) can be bounded 
from below by 

-C 5 N 2 F(A(N;p c ,p))F(A(N;p c , Pc )). 

This follows from a priori estimates of probabilities of two arms in a half- 
plane. We refer the reader to the proof of [20], Lemma 6.2, for more details. 
Again, although the proof of [20], Lemma 6.2, is given for site percolation 
on the triangular lattice, it also applies to bond percolation on the square 
lattice. 
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Putting together the bounds for the sums (40) and (41), and using (34), 
we obtain that the right-hand side of (39) is bounded from below by 

-C 6 N 2 F(A(N; Pc , P ))¥(A(N;p c , Pc )). 

Therefore, 

(44) ^log¥(A(N; Pc ,p)) > -CqN 2 F(A(N; Pc , Pc )) 

dp 

and 

F{A{N;p c ,p)) > P(i(Af;p c ,p c ))e- C6(p - pc)7v2lp( ^ (iV;Pc ' ? ' c)) 
>F{A(N;p c ,p c ))e- c ^ PN - p ^ N2¥ ^ N ' p ^ 
>C 7 F(A(N; Pc , Pc )). 
In the last inequality, we use (11). □ 

Definition 6.4. For any positive integers n <m < 2m < N and edge 
e £ Ann(m, 2m), we define C e (n,N;m) as the event that: 

- there exist two disjoint p c -open paths Pi and P2 inside Ann(n, N) \ {e}, 
the path Pi connects e x or e y to B{n) and the path P2 connects the other 
end of e to dB(N); and 

- there exists a p m -closed path P connecting e* and e* inside Ann(m, 2m)* \ 
{e*} so that PU {e*} is a circuit around the origin in Ann(m,2m)*. 

Note that if event C e (n, N;m) n {r e 6 (p c ,Pm)} occurs, then there is no p c - 
open crossing of Ann(n, N) and no p m -closed circuit in Ann(m,2m)* (see 
Figure 3). 

Definition 6.5. Let n, m and N be positive integers such that 2n < 
m and 3m < N. Let x = ([m/2], [3m/2]). For e G B(x,[m/2]), we define 
C e (n, N;m) as the event that: 

- the event A e (x;m;p c ,p m ) occurs; 

- there are two disjoint p c -open paths P5 and Pq such that P5 connects 
U\ m /2] (x) to the boundary of B(N) inside Ann(2m — 1, N) and Pq connects 
D[ m /2](x) to the boundary of B{n) inside Ann(n,m). Moreover, P5 and 
Pq satisfy P5 (~lAnn(x; [m/2], m) C (7[ m /2], m ( x ) and PenAnn(x; [m/2],m) C 

D[ m /2],m( x )'i 

- there exists ap m -closed path P inside Ann(m,2m— 1)* \B(x, [m/2])* such 
that P connects L[ m / 2 ](x)* to P[ m / 2 ](a;)* and P n Ann(x; [m/2], m)* C 

L[ m /2],m( x )* U ^[m/2],m( 2; )*- 

The event C e (n, N;m) is illustrated in Figure 4. 
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The event C e (n,N;m) obviously implies the event C e (n,N;m). The rea- 
son we introduce the event C e (n, N;m) is that 

(45) F(C e (n,N;m))~F(A e (x;m; Pc ,p m ))F cr (B(n)^dB(N)), 

where the constants do not depend on e, m, n and N . This observation 
follows from Lemma 6.1, the RSW theorem, and (35) and (36) applied to 
qi = Q2 = Pc- 

Corollary 6.2. For any positive integers n, m and N such that 2n < 
m and 3m < N , 

P(3e e Ann(m,2m):r e G (p c ,p m ),C e (n,N;m)) 

(46) 

>C s F cr (B{n)~dB(N)), 
where Cg does not depend on n, N and m. 

PROOF. Note that the events 

{r e G (pc,Pm),C e (n, N;m)} eeAnn ( m! 2 m ) 




B(n) 



Ann(m,2m) 



B(N) 



Fig. 3. Event C e (n, N; m) . The solid curves represent p c -open paths and the dotted curves 
represent p m -closed paths. The edge e does not have to be p c -open or p m -closed. 
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are disjoint. Therefore, 

P(3e € Ann(m,2m) :r e G (p c ,p m ),C e (n, N;m)) 
J2 F(T e (£(p c ,p m ),C e (n,N;m)) 

egAnn(m,2m) 

>(p m -Pc) J2 HC e (n,N;m)) 

e£B(x,[m/2]) 

>C 9 (p m - Pc ) Yl nMx;m-p c ,p rn ))¥ cr {B{n)^dB{N)) 
eeB(x,[m/2]) 

> C 10 ( Pm - p c )m 2 F(A(m;p c ,p c ))F cr (B(n) ^ dB(N)) 

>C n Wcr{B(n)^dB(N)). 

The second inequality follows from (45). In the third inequality, we use (37) 
and Lemma 6.3. In the last inequality, we use (11). □ 

Proof of Theorem 1.5. Lower bound. We give the proof for k = 2. 
The case k = 1 was considered in [18] and the proof for k > 3 is similar to 
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B(n)' 
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Ann(m,2m) 













B(N) 



Fig. 4. Event C e (n, N; m) . The solid curves represent p c -open paths and the dotted curves 
represent p m -closed paths. The edge e does not have to be p c -open or p m -closed. 
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the one for k = 2. Note that the event {P 2 > n} is implied by the event that 
there exists an edge e G B(n) and p> p c such that: 

- r e G (pc,p); 

- there exist two p c -open paths Pi and P 2 in P(n), the path Pi connects an 
end of e to the origin and the path P2 connects the other end of e to the 
boundary of B(n); 

- there exists ap-closed path P in B(n)* connecting e* to e* so that PU{e*} 
is a circuit around the origin. 

There could be at most one edge e G B(n) which satisfies the above three 
conditions. Therefore, 

[logn]-l 

P(P 2 >n)> F(3eGAnn([n/2 fe+1 ],[n/2 fc ]):r e G(p c ,p n/2fc+ i), 
fc=0 

C e (l,n;[n/2 k+1 ])) 

[logn]-l 

>Ci 2 ^ P cr (0^5P(n)) 
fc=o 

= Ci 2 [logn]P cr (0 «-» 95(n)). 
The last inequality follows from (46). □ 

7. Proof of Theorem 1.8. Let G = (G,£) be an infinite connected sub- 
graph of (Z 2 ,E 2 ) which contains the origin. We call an edge e G £ a discon- 
necting edge for G if the graph (Q, £ \ {e}) has a finite component and if the 
origin belongs to this finite component. Note that each outlet of the invasion 
is a disconnecting edge for the IPC. 

Let -D m>ri be the event that the IIC does not contain a disconnecting 
edge in the annulus Ann(m,n) and let T> m ,n be the event that the IPC does 
not contain a disconnecting edge in the annulus Ann(m,n). We prove the 
following theorem: 

Theorem 7.1. There exists a sequence (n^) such that 

(47) p(^I(P nfc , nfc+1 )<oo)=l 

and 

(48) u(Y l KD nkjnk+1 ) = o )=l. 

v k 1 
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Theorem 1.8 immediately follows from Theorem 7.1. Indeed, Theorem 7.1 
implies that the IIC is supported on clusters for which infinitely many of the 
events D nktTlk+1 occur and the IPC is supported on clusters for which only 
finitely many of the events T^n k ,n k+1 occur. Roughly speaking, this says that 
the distance between consecutive disconnecting edges (ordered by distance 
from the origin) can be much larger in the IIC than in the IPC. The proof 
of Theorem 7.1 is based on the following result (see Section 1.2 for the 
definitions). 



Theorem 7.2. There exist C\,Ci such that for all 1 <m < n, 

(49) P(£W) < CiPcrC^J 
and 

(50) u(D m>n )>C 2 ¥ ' 



■F cr (A- 



Lemma 7.1 ([17], Theorem 27). For all positive integers m <n and for 
all p G \pc,Pn], 

where the constants do not depend on m, n and p. 



Although Theorem 27 in [17] is stated for site percolation on the triangular 
lattice, the proof for bond percolation on the square lattice is the same. 



Lemma 7.2. There exists C3 such that for all m\ < ni2 < n, we have 

Proof. This follows from a priori estimates of probabilities of two arms 
in a half-plane (see [17], Theorem 24). □ 

Proof of Theorem 7.2. We first prove (49). Note that if the invasion 
percolation cluster contains a circuit, then there is a pond that entirely 
contains this circuit. Therefore, the event T> m ^ n can only occur if there exists 
an invasion pond which contains two disjoint crossings Pi and P2 of the 
annulus Ann(m,n) (see Figure 5). Therefore, there exists p' such that P\ 
and P2 are p'-open and there exists a circuit around the origin which is 
p'-closed and which has diameter at least n. 
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Ann(m,n) 



Fig. 5. Event X> m ,n. The edges e and f are disconnecting. The paths Pi and P2 create 
a circuit, which implies that there is a pond that entirely contains both paths. 



Recall the definition of (p n (j)) from (5). Later, we take C* in (5) to be 
sufficiently large. We decompose the event T> m n according to the value of 
p<: 

log* n 

(51) P(X»m,n) = £ HVm,n,p' € \p n (j),Pn(j ~ I)))- 

5=1 

Note that the event {T> m>n ;p' G [p n (j),Pn(j — 1))} implies the event 
A mn p„(j-i) ^ B n ,p n (j) ( see Section 1.2 for the definition of these events). 
It follows from (6) and (10) that there exist constants C4 and C5 such 
that the probability P(i? n pn (j)) is bounded from above by C^log^ -1 - 1 n)~ C5 . 
The constant C5 can be made arbitrarily large by making C* large enough. 
We use Lemmas 7.1 and 7.2 to bound the probability 1P(^4^ n Pn (j-\)) ^ 

C6(log^ -1 ^ n)F cr (Am n ). We use the FKG inequality and the above esti- 
mates for the events ^m,n, Pn (i-i) and B n, Pn (j) to S et 

log* n 

nV m , n )<C 4 C 6 F cr (A 2 ^ n ) £ (log^M^ 5 

i=i 

< C 7 ¥ cr (A^ n ). 

The last inequality follows from [6], (2.26), if we take C* such that C5 > 1. 
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We now prove (50). Let C m , n be the event A[ m / 2 }, m n A \ m /2],2n n A ™M (see 
Figure 6). Note that v{C m , n \ f m , n ) = 0. It is therefore sufficient to prove 
(50) for C m 

,71 ■ 

For positive integers m < n < N (later, we consider the limit as iV tends 
to infinity), we use the FKG inequality to get 

> P cr (C min , n A 0rn n A^ N ) 

> C8P > cr(-4f m /2],2n) P cr(-4o, m )I !) cr(^n,Ar) 

(A 



>c 8 



~^cr(A mn ) 



for some Cg > 0. Standard RSW arguments give a constant Cg such that for 
all 1 < m < n, 



cr(^f m /2],2n) — C^cr^m^)- 

Therefore, 





v{D ^ n) - c ^^Afty □ 
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Proposition 7.1. There exists a sequence (rife) such that rik+i >An k , 

( 52 ) E P -«,n fc+1 )<^ 

k 

and 



(53) E 



OO. 



Proposition 7.1 follows from Lemma 7.2 and the fact that F cr (Al nn ) < 
c(m/n) s for some positive c and 8. Indeed, we obtain W cr (A]. nn ) < 
c(C^S' cr {A^ m n )Y ■ We now take, for example, the sequence n k = min{ra > 
4n fc _i:P cr (A2 fc _ ijn )<(l/fc)i+*}. 

Proof of Theorem 7.1. We take a sequence from Proposition 7.1. 
Equality (47) follows from the Borel-Cantelli lemma. To prove (48), we use 
Borel's lemma [13]: 

Lemma 7.3. Consider a probability space {Sl,J~, P) and a sequence of 
events T n € T . Let limsup n r n = C\n^Jk>n^k be the event that infinitely 
many of the T n 's occur. Let a n =I(T n ) be the indicator of event T n . If 
there exists a sequence b n such that ^l n b n = oo and for any en € {0,1}, 
i = 1, . . . , n — 1, 

P(r n |ai =ai,...,o n _ 1 = a n _i) >b n > 0, 

then 

Pllim sup T T: 



Note that it is sufficient to prove (48) for the events C nkt n k+1 (see the proof 
of Theorem 7.2 for the definition). We apply Lemma 7.3 to the probability 
measure v and to the events C n2k ^ 2k+1 . Let d k = l{C n2k ^ 2k+1 ). A slight 
extension of the proof of (50) gives, for any on £ {0, 1}, i = l,...,k — l, 

(54) V{C n2k ,n 2k+l \dl =«!,..., 4-1 = OL k -i) > C 2 7~A1 \ =: 

where C2 is the constant from (50). Indeed, let W be the set of configurations 
of edges in B{2ri2k-i) such that d\ = a±, . . . ,4-1 = a k-i- F° r an Y w G W 
and large enough N, 

~^cr{Cn 2k ,n 2k+1 D A\ N \uj) 

> ^cr(C n2k . n2k+1 n Al n2k n A\ 2 n \uj) 
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"™2fc+i 



,n\u) 



which implies (54). In the second line, we used the FKG inequality and 
independence. The equality follows from independence. Prom the choice of 
(nfc), it follows that J2k — °o. Therefore, equality (48) follows from Lemma 
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